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Abstract. We investigate pattern speeds in spiral galaxies where the structure is induced by 
an interaction with a companion galaxy. We perform calculations modeling the response of 
the stellar and/or gaseous components of a disc. Generally we do not find a unique pattern 
speed in these simulations, rather the pattern speed decreases with radius, and the pattern 
speed for individual spiral arms diff'er. The maximum pattern speed is ~20 km s~^ kpc~^ for 
the discs with a live stellar component, decreasing to 5 km s~^ kpc~^ at the edge of the spiral 
perturbation. When only the gas is modeled, Qp is typically very low (5 km s~^ kpc~^) at all 
radii. 
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1. Introduction 

Most grand design spiral patterns are believed 
to b e due to interactions wi th other galax- 
ies (iToo mre & Toomrd 1972h. or are driven 
by ba rs ( Korme ndy & Normanll 19791: iBottemal 
12003, Athanassoula et al., this volume). Here 
we investigate the first of these scenarios, by 
performing numerical simulations of interact- 
ing galaxies. In this instance, there is not nec- 
essarily a singular pattern speed for the spi- 
ral arm s, as indicated by recent observations of 
M51 (Meidt et al.ll2QQ8ah . 

2. Method 

We use the Smoothed Particle Hydrodynamics 
Code (SPH) to model a galaxy subject to an 
interaction. The first galaxy is modeled by as- 
suming a spherical potential for the halo, with 
3 diflferent scenarios for the galactic disc: a) a 



live stellar disc, no gas; b) a logarithmic stel- 
lar potential and a live gaseous disc and c) a 
disc containing live gaseous and stellar com- 
ponents. The stars and/or gas are allocated ve- 
locities such that the Toomre instability param- 
eter, Q, is globally 2, but in addition, the galaxy 
was allowed to evolve in isolation until any 
flocculent structure disappears. In all simula- 
tions 1 million particles are used. For the case 
of stars and gas there are 500,000 gas particles 
and 500,000 stellar particles, and in all cases 
the total mass of the disc is 5 x 10^ M©. 
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We adopt a similar approach to lOh et al.l 
(l2008h to model the interaction. The inter- 
acting galax y is repr esented by a sink parti- 
cle ( Bate et. a l] ll995h and is of relatively low 
mass, equal to the mass of the disc or 2.5% of 
the total mass of the first galaxy. The interact- 
ing galaxy takes a parabolic orbit, reaching a 
closest approach of 25 kpc after a time of 370 
Myr. Initially this galaxy is at a distance of 50 
kpc. 
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Fig. 1. The column density is shown for a stellar 
disc of a galaxy perturbed by a (point mass) galaxy 
on a parabolic orbit. The second galaxy is located 
70 kpc away, at the coordinates (-20, 68) kpc. The 
stellar disc contains broad and diffuse spiral arms. 



Fig. 3. The column density is shown for the gas disc 
perturbed by a galaxy, at the same time as Fig. 1. 
The gas shocks produce very narrow, dense spiral 
arms compared to Fig. 1. 
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Fig. 2. The pattern speeds for each spiral arm are 
plotted versus radius for the stellar disc, is dif- 
ferent for each arm and decreases with radius. The 
lower line shows a 1 /r dependence. 



25 : 

20 - 



- Arm 1 

- Arm 2 



5 10 
Radius (kpc) 



Fig. 4. The pattern speeds are plotted for the spiral 
arms of the gas disc. Qp is much lower than typical 
from observations. 



3. Results 

3.1. Stellar disc 

Fig. 1 shows the column density of a stellar 
disc, at a time of 820 Myr. The interacting 
galaxy is 70 kpc from the centre of the plot. 
The spiral arms are relatively weak and broad. 

We calculate the pattern speed of the spiral 
arms according to 

^ _ 0(Pmax)l?2 - 0(Pmax)k ... 



where pmax is the peak density of a particu- 
lar spiral arm at a given radius, and ti and t2 
are times during the simulation. First we se- 
lect points covering a particular spiral arm at 
time ti = 800 Myr, when a strong spiral pat- 
tern has emerged. Then the azimuthal angle of 
the spiral arm is located for different radii. This 
process is repeated at time t2 = 925 Myr, to ob- 
tain the change in azimuthal angle of the spiral 
arm at each radius, and thus Qp. Given there 
are 2 spiral arms, this method leads to a pat- 
tern speed for each spiral arm. These pattern 
speeds are shown versus radius in Fig. 2, the 
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errors reflecting the uncertainty in locating the 
peak density of the spiral arm. 

The magnitude of the pattern speeds are 
not dissimilar from som e of those measured 
for spirals (e. g., Clemens "& Alexandeill2QQll : 
lGrosb0l et al.| [2bQ6). although the spiral arms 
clearly exhibit difl'erent pattern speeds. This 
difl'erence is a consequence of the asymmetry 
of the system, i.e. that the interaction induces 
one arm on one side of the disc first. The pat- 
tern speeds for each arm also decrease with ra- 
dius, roughly as expected for spiral patterns in- 
duced by interactions. 

3.2. Gaseous disc 

We also performed calculations with just gas. 
The gas constitutes 1% of the mass of the 
galactic disc. This value is unrealistically 
small, but the low gas mass is chosen to avoid 
gravitational instabilities, which would halt the 
calculation. Essentially, we are only investigat- 
ing the reaction of the gas to the interaction, not 
the self gravity of the gas. These calculations 
are also not a particularly realistic case as they 
ignore the perturbation experienced by the stel- 
lar disc from the interaction (instead the stellar 
disc is represented by a symmetric potential). 
However the case with just gas is explored for 
completeness. 

In Fig. 3, we show the disc when only the 
gas is included, at the same time (800 Myr) as 
Fig. 1. The spiral arms are clearly much nar- 
rower, more prominent, and more dense than 
for the stellar disc. The evolution of the gas and 
stellar discs is also diff'erent. The gaseous spiral 
arms rotate much slower than the stellar arms. 
Consequently the pattern speed is very low, ~4 
km s"^ kpc"^ and does not show much varia- 
tion with radius (Fig. 4). The diff'erence com- 
pared to the stellar disc is that one spiral arm is 
still linked to, and rotates at the same angular 
velocity as, the orbiting galaxy. 

3.3. Evolving botli the stars and gas 

Finally we show calculations with gas and 
stars, again where the gas represents 1% of 
the mass of the stellar disc. Thus the mass of 



Fig. 5. The gas distribution is shown from a cal- 
culation with gas and stars at times of 500 (top), 
740 (middle) and 1060 Myr (bottom). The spiral 
arms are very asymmetric at the earliest time (when 
the position of the interacting galaxy is located at 
(22.5, 15) kpc), but becomes more symmetric at 
later times, after the interaction. 



the stellar disc is 5 x 10^ M© and the gaseous 
disc 5 X 10^ Mq, although the actual num- 
ber of gas and stellar particles in the calcula- 
tion are equal. Fig. 5 shows the column den- 



110 



Dobbs: Interacting galaxies 




- Arm 1 

- Arm 2 



5 10 15 

Radius (kpc) 

Fig. 6. The pattern speed is shown for the gaseous 
arms, from the simulation with stars and gas. The 
gas is largely coupled to the stars, hence Qp is higher 
than when only gas particles are present. Again there 
is a clear difference in the pattern speeds of each 
spiral arm. 



sity of gas from the calculation with live stellar 
and gaseous components. The morphology evi- 
dently changes over time. At an earlier stage in 
the interaction, the spiral pattern is very asym- 
metric (unlike the calculations with gas alone, 
where the spiral pattern is symmetric through- 
out the simulation). The pattern becomes more 
symmetric at later times, tending towards the 
distribution without a live stellar disc (Fig. 3). 
Although not shown on Fig. 5, the stellar distri- 
bution reflects the gas distribution, but the spi- 
ral arms are much weaker and broader. 

Finally in Fig. 6 we plot the pattern speed 
of the gas from the calculation with both live 
stellar and gaseous components. The times se- 
lected to calculate the pattern speed are 450 
and 500 Myr. The spiral arms clearly have very 
difl'erent pattern speeds, which is not surprising 
given the asymmetry of the disc. The spiral arm 
with the lower pattern speed is still associated 
with the interacting galaxy. The pattern speeds 
are more similar in magnitude to the case when 
only stars are used (Fig. 2), since the gas is bet- 
ter coupled to the stars. 

4. Conclusions 

We have performed calculations with a stellar, 
gaseous and both stellar and gaseous disc sub- 
ject to an interaction with an orbiting galaxy. 



The pattern speed across the disc is generally 
not constant, and pattern speeds in each arm 
difl'er. For a stellar disc, Qp = 5 - 20 km s"^ 
kpc"\ decreasing with radius approximately as 
1/r. With only gas, the pattern speeds are much 
lower (3-6 km s"^ kpc"^). When stars and gas 
are included, the gas tends to follow the stel- 
lar distribution, thus the pattern speeds of the 
gaseous spiral arms are higher (5-17 km s"^ 
kpc-i). 

These calculations may be improved by 
using a more consistent initial g alaxy set up 
(e.g., iKuiiken & Dubinskil Il995l) . A natural 
extension of this work would also be to 
compare with observa tions by applying the 
Tremaine- Weinberg (iTremaine & WeinbergI 
1984) method to these calculations (see also 
Meidt et al. 2008b). 
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